Disentangling non-Gaussianity, bias and GR effects in the galaxy distribution 
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Local non-Gaussianity, parametrized by /nl, introduces a scale-dependent bias that is strongest at 
large scales, precisely where General Relativistic (GR) effects also become significant. With future 
data, it should be possible to constrain /nl = C(f) with high redshift surveys. GR corrections to 
the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to 
/nl = O(l), so it is essential to disentangle these effects. For the first time in studies of primordial 
non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the 
importance of a proper definition of bias. We present observable power spectra with and without 
GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, 
these effects can be distinguished by their different redshift and scale dependence, so as to extract 
the true primordial non-Gaussianity. 



I. INTRODUCTION 

The galaxy power spectrum has proven an indispens- 
able tool in understanding our Universe, and measur- 
ing it underlies the motivation for many future surveys. 
Usually the 3D power spectrum P{k) is inferred from 
the measured galaxy redshifts and angular positions and 
compared to predictions based on constant-time slices 
of the Universe. However, surveys are slices down the 
past lightcone of the observer; constant-time slices are in 
principle unobservable for significant redshifts and, for a 
given observer, are not uniquely associated with redshift. 
Furthermore, for large-scale perturbations of the order of 
the Hubble radius H , P(k) needs to be defined for the 
matter overdensity S m in a given choice of time slice, i.e. a 
choice of coordinates or "gauge" , which can be arbitrarily 
specified. On sub-Hubble scales, the power spectrum is 
not sensitive to the choice of gauge, but significant am- 
biguities arise near and beyond H^ 1 . It is possible to 
define a gauge- invariant 5 m , but this does not solve the 
problem. Instead, one must compare predictions for true 
observables, such as angular power spectra at a given 
redshift, which are necessarily gauge invariant. 

In order to relate the observed overdensity to the over- 
density defined in any chosen gauge, General Relativis- 
tic (GR) corrections must be included to take account 
of both gauge effects and lightcone effects. These correc- 
tions have recently been analyzed in detail [H-Q ■ The cor- 
rections are relevant close to and beyond the Hubble scale 
H -1 . For highcr-rcdshift surveys, the corrections also be- 
come significant on smaller scales, since H~ 1 (z) decreases 
with the redshift z. For example, H~ 1 (z = 2) ~ QAHZ , 
and H~ x (z = 10) ~ O.ObH^ 1 . The angular extent of the 
Hubble radius at redshift z is 9h{z) = H~ 1 (z)/Da(z), 
where Da(z) is the angular diameter distance. At z = 2, 
we have 6h ~ 50°. 

Just as the overdensity definitions depend on the 
gauge, so too does the galaxy bias, which relates the fluc- 
tuations of galaxy number density to the underlying mat- 
ter density fluctuation. An ambiguity in the definition of 
bias can not only confuse our understanding of the un- 



derlying density power spectrum, but it also can impact 
our understanding of the non-Gaussianity of the under- 
lying density field, because local type non-Gaussianitics 
characterised by /nl have been shown to produce a scale- 
dependent bias @, II] (for a review see Ref. [j| and refer- 
ences therein). 

On large scales, measurements are limited by cosmic 
variance. However, a recent proposal aims to eliminate 
cosmic variance by correlating a highly biased tracer of 
large-scale structure against an unbiased tracer 0. Us- 
ing this method, constraints at the level /nl = O(l) are 
estimated [HI- The GR corrections and the gauge de- 
pendence of the bias introduce effects in the observed 
angular power spectra of the same order as the scale de- 
pendent bias effects with /nl = 0(1), as we show below. 
Thus it is essential to include the GR corrections and 
clarify the gauge dependence of the bias prescriptions in 
order to obtain optimal constraints on primordial non- 
Gaussianity. 



II. PRESCRIPTIONS FOR GALAXY BIAS 

It is generally assumed that the galaxy bias is scale- 
independent on large scales if the primordial fluctuations 
are Gaussian. However, once GR corrections become 
important, it is essential to specify the physical frame 
in which the scale-independent bias is defined. In [D-Q 
the scale-independent bias is defined in the uniform red- 
shift gauge, i.e. it relates the observed galaxy overdensity 
to the matter density in the same observational gauge. 
However, in general bias is determined by local physics 
that governs the formation of galaxies and is inde- 
pendent of how they are observed. 

Local bias should be defined in the rest frame of 
CDM (which is typically assumed to coincide with the 
baryon rest frame). Bias can be understood in the peak- 
background-split approach where halo collapse occurs 
due to small scale peaks in the density exceeding a critical 
value, 5+ [13j. Large-scale fluctuations <5i ong lead to an en- 
hancement of the density of collapsed halos by effectively 
reducing the critical value for small-scale fluctuations to 
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— S\ ong . In the spherical halo collapse model, there is 
an exact GR interpretation of the standard Newtonian 
collapse [3]. The criterion for collapse of a local over- 
density is that the linear overdensity in the comoving- 
synchronous gauge reaches the critical value, <5* = 1.686. 
Thus linear scale- independent bias applies on large scales 
for a Gaussian distribution in the comoving-synchronous 
gauge (as used by |4|-l6j|). By contrast, in the Newtonian 
(longitudinal) gauge there is no local criterion for collapse 
in terms of the overdensity on large scales 14j. 

In comoving-synchronous gauge (ignoring vector and 
tensor modes) 

ds 2 = a(7 ? ) 2 {-<V+ {l-2K)5 ij + 2d i d j E t&W'}, (1) 

and the CDM velocity perturbation vanishes in this 
gauge. In the simple linear model, the Eulerian bias b 
is defined via 



where Sn 
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Sn^/n and 8m = Sp^/p are the 
galaxy and matter fractional overdensities in comoving- 
synchronous gauge. 

Although the observed galaxy overdensity can be cal- 
culated in any gauge, it is often convenient to use the 
Newtonian gauge because it facilitates comparisons with 
the standard Newtonian calculations. In Newtonian 
gauge, 



ds 2 = a{r]) 2 - (1 + 2<J>)cfy 2 + (1 - 2^)dx 2 



(3) 



In the matter-dominated era, we can set \t = $ since the 
effect of radiation anisotropic stress is negligible. From 
now on, quantities without superscripts refer to New- 
tonian gauge values if there is ambiguity. The density 
perturbations in the two gauges are related by 



n k 



p K 



where v is the velocity. Thus in Newtonian gauge, using 
@ we have S n = bS m — [b(\np)' — (Inn)'] v/k. Assuming 
that galaxy number and matter are both conserved, we 
have n'/n = p'/p = —ZaH and hence 
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(5) 



On small scales, a/k <C H^ 1 , the gauge difference van- 
ishes and S ~ 5( c > ~ <5 std , where ^ std is the standard 
density perturbation in the Newtonian limit. Note that 
for the growing mode in the matter era (i.e. neglecting 
any dark energy effect), in Newtonian gauge we have 
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where the Newtonian potential $ remains constant on all 
scales. In comoving-synchronous gauge 
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The two gauges coincide on sub-Hubble scales, but they 
differ significantly near the Hubble scale. 

In 0Ta|; t ne analysis of GR corrections to the power 
spectrum is performed for a Gaussian primordial spec- 
trum. Here we generalize to include primordial non- 
Gaussianities of the local form. The potential $ has a 
non-Gaussian contribution: 



$ = ^-/nl(^-(^)), 



(8) 



where (j) g is Gaussian 



The Newtonian potential $ is 
related to the density perturbation in the comoving- 
synchronous gauge by the Poisson equation ([7]). There- 
fore this local form for the Newtonian potential corre- 
sponds to a non-local form for the density field in the 
comoving-synchronous gauge. Following [7, 8], this non- 
Gaussianity induces a scale-dependent correction to the 
bias: 



6 + /nl0-1) 



k 2 T(k)D(a) 



(9) 



where b is the Gaussian bias, T(k) is the transfer function 
(fc-independent at low k), and D(a) is the growing mode 
of density perturbations. We emphasize that this defi- 
nition is intrinsically relating the comoving-synchronous 
gauge densities [14|. 

In this paper, we only include local-type primordial 
non-Gaussianity ([S]) for simplicity. However, any non- 
Gaussianity that produces a non-zero bispectrum of the 
Newton potential in the squeezed limit in the matter era 
leads to a scale dependent correction to the bias (l5j . 
This includes non-linear evolutions of the potential on 
super- horizon scales [16( and it should be added to ([9|). 



III. GAUGE DEPENDENCE OF THE POWER 
SPECTRUM 

FigureQ] (left panel) illustrates how the number density 
power spectrum P{k) depends on the gauge choice of the 
variables, on the gauge used to define the bias and on the 
non-Gaussianity. The simple linear behaviour (solid line) 
results from plotting the comoving-synchronous gauge 
number density using a linear scale-independent bias de- 
fined in comoving-synchronous gauge, ([2]). The corre- 
sponding Newtonian gauge number density (dashed line) 
differs by a term 0{k~ 2 ). This scale dependence is ex- 
actly the same as the scale dependence of the bias in 
the comoving-synchronous gauge arising from primordial 
non-Gaussianity ©. 

Consider the galaxy overdensity 5 n given by ([5]) in the 
matter era. On large scales, the velocity term 3aHv/k 
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FIG. 1: Left: The power spectrum of various density perturbations at z = 1. We use a standard LCDM background and 
assume b — 2. Right: The angular power spectrum at z — 1 assuming all galaxies are in a Gaussian window function of width 
o z = 0.1. See text for detailed explanation. Note that we take /nl = 0.63, which is slightly different from /nl = 3/^4 due to 
the effect of the cosmological constant. Also A evolves with redshift and a different /nl would be required at another redshift. 



and the non-Gaussian correction to the bias © dominate 
and S n behaves as 



3 + /nU(6 - 1) 



( Qg ) 5^ (10) 
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where we defined 
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where $i n it is the potential in the radiation era. The 
gauge correction to the galaxy overdensity in Newtonian 
gauge, corresponding to the first term in square brack- 
ets in (ITU1) . gives a similar 0(k~ 2 ) effect as the non- 
Gaussian correction to the bias, the second term. In 
the matter-dominated era, the power spectrum of the 
Newtonian gauge 5 n coincides with that in the comoving- 

synchronous gauge 5n \ with /nl = 3/[A(b — 1)]. We em- 
phasize that this is an apparent effect - if we calculate 
true observables, they are independent of the choice of 
gauge and there is no confusion. However, this clearly 
demonstrates the importance of using true observables 
that include all the GR corrections when analyzing the 
effect of primordial non-Gaussianity. 

In addition to the gauge dependence of the galaxy 
overdensity, the ambiguities in the gauge used to define 
bias introduce another potential confusion to the scale- 
dependent bias effect from primordial non-Gaussianity. 
The dotted line shows the impact on the Newtonian 
gauge overdensity of wrongly defining a linear bias in 
Newtonian gauge, 5^ Tone — bS m . On large scales, this 
behaves as 



jwrong ^ 3^ 



(aHY 
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5< c >. 



(12) 



Comparing (fT0|) and (fT2|) . we find that if the linear bias 
is assumed to be scale-independent in Newtonian gauge, 
this would be interpreted as the effect of primordial non- 
Gaussianity with /nl = 3/ 'A in the matter-dominated 
era. 



IV. OBSERVED ANGULAR POWER SPECTRA 

Figure [T] (left panel) illustrates how non-Gaussianity 
could be confused in the power spectrum with GR gauge 
corrections on large scales. However, the power spectrum 
and the overdensities themselves are not directly observ- 
able. In the Newtonian limit, it is standard to define the 
observed galaxy overdensity as 



sld 



dv 



(13) 



where <5^ td is the number density fluctuations in the New- 
tonian limit, v the peculiar velocity, and k is the lensing 
convergence 0] 



07 - v) 



v. r 

2 n J Vo (Vo-v)(Vo-fj) 



-($ + V)dfj. (14) 



However, there is an ambiguity in the definition of 5„ in 
this formula on large scales because we need to specify the 
gauge for 5^ d . We must also include all effects impacting 
the observed galaxy overdensity, including redshift-space 
distortions, gravitational lensing and integrated Sachs- 
Wolfe effects. 

The observed galaxy overdensity including these ef- 
fects has been first derived in [J Q. In these papers, 
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bias is defined on the constant observed redshift hyper- 
surface. As discussed in section III, bias should be de- 
fined in the comoving-synchronous gauge and a wrong 
definition of bias introduces spurious scale dependent ef- 
fects in the power spectrum. Refs. [5j, [6| derived the ob- 
served galaxy overdensity in Newtonian gauge by cor- 
rectly defining bias in the comoving-synchronous gauge. 
The observed galaxy overdensity in the direction h at 
redshift z is obtained as 



A n (n,z) = S n 



H 



1 » dv ^ 
aH d\ 



H 2 aH X 



X 



<7o l 

($ + ^)dr] + — + $ - 2* 
aH 



J ($' + *')d?7 - n v 
(15) 



where all perturbations are expressed in Newtonian 
gauge, a subscript 'o' denotes the observer, and x is the 
comoving radial position of the source. The first three 
terms are the same as the standard formula in the New- 
tonian limit. The second line comes from the redshift per- 
turbation of the volume, which contains a Doppler term, 
the ordinary and integrated Sachs- Wolfe terms. The in- 
tegral in the third line comes from the radial shift due 
to lensing. The potentials, ^ and are suppressed on 
small scales compared with density fluctuations but they 
become the same order near the horizon scales and their 
effects need to be taken into account consistently. 

In order to compute the observable effects of GR cor- 
rections and non-Gaussianity, we consider the angular 
power spectrum Ci of A n in a redshift slice. The counts 
in a redshift slice are given by 

A n (h,z) = y^ j at m {z)Y im (h), C e (z) = {\a em (z)\ 2 }. 

This angular power spectrum includes the auto- and 
cross-correlations amongst all terms on the right of (fl~5|) , 
taking into account the relevant transfer functions 

Figure Q] (right panel) shows the angular power spec- 
trum at z = 1 assuming all galaxies are in a Gaussian 
window function of width a z =0.1. As emphasized be- 
fore, the observed angular power spectrum is indepen- 
dent of gauge, so that there is no ambiguity in the choice 
of gauge for the galaxy density perturbation. However, 
there still remain the ambiguities in the gauge used to 
define bias. The dotted line shows the result of defining 
a linear bias in the Newtonian gauge, which is degenerate 
at fixed redshift with a scale-dependent bias effect from 
primordial non-Gaussianity. 

In Figure [3J we compare the full angular spectra with 
the standard formula based on the Newtonian limit re- 
sults, Eq. (1131) . As mentioned above, there is an ambigu- 
ity in the definition of <5„ in this formula on large scales. 
We use the definition <5^ td = bd$ because the growth rate 

(c) 

of 6m is the same as the Newtonian result on all scales 



and the gauge difference is suppressed by the 1/fc factor 
in Eq. © 

1- 

The top-left plot shows the impact of full GR correc- 
tions on observed spectra compared to the standard re- 
sult at various redshifts, for Gaussian fluctuations. At 
higher redshift, the GR effect becomes larger and more 
important at larger £, since the Hubble scale is smaller. 
These corrections are suggestive of the changes which 
arise due to primordial non-Gaussianity, which are shown 
in the top-right panel, including all GR corrections. How- 
ever, at high redshifts, we find that the shape of the angu- 
lar power spectra induced by the GR corrections is very 
different from the shape induced by the non-Gaussian ef- 
fect, and these two effects look very different at small 
I. This is demonstrated in the bottom-left panel, where 
we attempted to tune /nl in the standard formula so 
that it cancels GR corrections. This lack of degener- 
acy implies that it is in principle possible to disentangle 
between the two effects in the observed angular power 
spectrum at different redshifts. The bottom-right panel 
shows the GR corrections in the presence of primordial 
non-Gaussianity /nl = 50. At low redshifts, the non- 
Gaussian effect is still small and GR corrections are of 
the same order as the Gaussian case. As expected, GR 
corrections become unimportant for small t and at high 
redshifts where the non-Gaussian effects dominate. (All 
numerical computations have been performed with the 
freely available package CMBquick jl7|].) 

Note that in our analysis we have considered the sim- 
plest case of conserved galaxy number and neglected se- 
lection effects, e.g., luminosity dependence in a realistic 
magnitude-limited survey @] . It is also interesting to 
consider possible degeneracies with other physical effects 
that have an influence on the power spectrum on large 
scales, such as clustering dark energy, tilted spectrum 
and isocurvature modes, which would need further inves- 
tigation. 



V. CONCLUSIONS 

We observe on our past lightcone and the correct inter- 
pretation of data on very large scales requires GR correc- 
tions to the standard Newtonian approach. This is par- 
ticularly important for detecting non-Gaussianity using 
large-scale structure. For the first time, we have included 
all GR corrections to the galaxy power spectrum in an 
analysis of primordial local-type non-Gaussianity. We 
only considered local-type primordial non-Gaussianity 
(|8|) in this paper because any other types of primordial 
non-Gaussianity that produce a negligible bispectrum of 
the Newton potential in the squeezed limit do not lead 
to a scale dependent correction [ll|. We have shown 
that galaxy bias is naturally defined in the comoving- 
synchronous gauge, and choosing a different gauge can 
lead to incorrect conclusions on scales approaching the 
Hubble radius where GR corrections become important. 



5 



1.02 



1.01 



1.00 



0.99 



0.98 



0.97 



1.01 



0.97 




loo _ -z-X— 



0.99 



0.98 



- 0.98 



FIG. 2: Ratio of the standard angular power spectra to the full GR spectrum in various cases. See text for detailed explanation. 
We assume 6 = 2 and all the galaxies are in a Gaussian window function of width <r z — O.lz. 



GR corrections to the power spectrum and ambigu- 
ities in the gauge used to define bias introduce effects 
qualitatively similar to local-type non-Gaussianity with 
/nl = However we have shown, using the ob- 

served angular power spectrum, that while the wrong 
bias definition (|12p has the same scale-dependence as 
non-Gaussian effects, it has different redshift depen- 
dence. Moreover, GR effects have a somewhat different 
scale dependence. These differences arc the basis for dis- 
criminating between these effects and identifying the true 
non-Gaussian signal. 

GR effects are important at wide angles and at high 
redshifts as these effects are important near horizon 
scales (see Fig. 1 and Fig. 2). Thus a deep and wide angle 
survey makes the effects easier to see. On large scales, 
measurements are limited by cosmic variance. Given 
that GR corrections are comparable to local-type non- 
Gaussianity with /nl = they are not measurable 
without reducing cosmic variance. However, a recent pro- 



posal aims to eliminate cosmic variance by correlating a 
highly biased tracer of large-scale structure against an 
unbiased tracer (Toj . It is important to study whether 
general relativistic effects are measurable in future sur- 
veys [H. 
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